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Abstract. This review starts with a discussion of the hierarchy of scales, relevant to 
the description of superfluids in neutron stars, which motivates a subsequent elemen- 
tary exposition of the Newtonian superfluid hydrodynamics. Starting from the Euler 
equations for a superfluid and a normal fluid we apply the tensor virial method to 
obtain the virial equations of the first, second, and third order and to compute their 
Eulerian perturbations. Special emphasis is put on the computation of perturbations of 
the new terms due to mutual gravitational attraction and mutual friction between the 
two fluids. The oscillation modes of superfluid Maclaurin spheroids are derived from 
the first and second order perturbed virial equations. We discuss two generic classes 
of oscillation modes which correspond to the co-moving and relative oscillations of two 
fluids. These modes decouple if the normal fluid is inviscid. We also discuss the mixing 
of these modes (when the normal fluid is viscous) and its effect on the dynamical and 
secular instabilities of the co-moving modes and their damping. 

1 Introduction 

Radio and x-ray observations of neutron stars provide strong evidence for the 
superfluidity of neutron star interiors. Perhaps, the most striking manifestations 
of their superfluidity are the long (on time-scales from several hours to hundreds 
of days) relaxations that follow the glitches in the spin and spin-down rates of 
some pulsars. Although the majority of pulsars are very precise clocks, timing 
observations reveal persistent random fluctuations in times of arrival of radio 
signals. Some pulsars show long-term periodicities in their spin-characteristics 
and periodic changes of their pulse shape. However, the relation between the 
superfluidity of neutron star interiors and these anomalies of pulsar timing is 
not firmly established yet. 

Further evidence for the superfluidity of neutron star interiors came in the 
1980s with the advent of the orbiting x-ray satellites. The measurements of 
the thermal radiation from a dozen or so hot neutron stars provided indirect 
information on the temperature of superfluid phases in neutron stars. Theoretical 
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thermal histories of superfluid neutron stars are consistent with the x-ray data 
(within the limits of our knowledge of the input physics). Non-superfluid stars, 
as a rule, cool too fast to below the threshold of detection. 

Apart from the radiation in the electro-magnetic spectrum, neutron stars 
are expected to be primary sources of gravitational wave radiation, which are 
expected to be detectable by future laser interferometer detectors. It is hoped 
that one can probe neutron star interiors and their superfluidity using gravity 
waves, as their eigen-frequencies and damping may depend on the dissipation in 
the superfluid, at temperatures below the superfluid phase transition. 

This review concentrates on the oscillations of superfluid self-gravitating el- 
lipsoids within the tensor virial method. The method was originally introduced 
by Chandrasekhar and Fermi in the context of magneto- hydrodynamics [jjj. It 
was extensively developed in the 1960s by Chandrasekhar for the study of the 
ellipsoidal figures of equilibrium and their oscillations. A comprehensive account 
of this work is contained in Chandrasekhar's monograph Ellipsoidal Figures of 
Equilibrium (hereafter EFE)^]. 

The ellipsoidal approximation provides an idealized picture of oscillations of 
neutron stars. One can think of several arguments in favor of adopting such an 
approach: first, the combination of the Newtonian gravity and two- fluid hydro- 
dynamics of superfluid defines an exactly solvable model if we assume that the 
fluids are incompressible and inviscid; second, past experience with single-fluid 
self-gravitating ellipsoids shows that most of the qualitative features found for 
these ellipsoids have their analogs in more "realistic" systems; third, the method 
is transparent and in many cases analytical results can be obtained which shed 
light on the underlying physics. 

The tensor virial method is not the only tool for investigating the properties 
of ellipsoidal figures. Alternative formulations exist in the literature and we refer 
to Refs. ||,(|||] for further details; for a pedagogical introduction see the textbook 
[pi. Note, also, that various formulations of the theory of ellipsoids, to a large 
extent, are equivalent. 

Superfluid oscillations were studied using various methods and approxima- 
tions in the past decade. Epstein pointed out that the superfluidity of neutron 
star interiors has potentially important effects on the propagation of seismi- 
cally excited acoustic waves. It allows for additional types of waves to propagate 
by virtue of doubling of degrees of freedom in a superfluid; superfluid phases 
create acoustic discontinuities in which wave velocities or polarizations change 
abruptly on the bounding interfaces JtJ. The effects of superfluidity on global 
hydrodynamic oscillations were investigated by Lindblom and Mendell S in a 
model where the superfluid and the normal fluid are coupled via gravitational 
attraction and the entrainment effec^. Their solutions reveal that the lowest 
frequency pulsations are almost indistinguishable from those derived from the 
ordinary-fluid hydrodynamics; however, their analytical solutions also reveal the 

1 The latter effect arises in the layers where neutrons and proton condensates coexist, 
hence the flow of one condensate is accompanied by the motion of the other. See 
B. Carter's article in this volume for further details. 
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existence of a spectrum of modes which are absent in a single fluid star. Non- 
radial oscillations of non-rotating superfluid neutron stars were computed by 
Lee, whose numerical solutions for the radial and non-radial pulsations of the 
two-fluid stars identified distinct superfluid modes ||. The effects of shear vis- 
cosity of the electron fluid and mutual friction on the r-mode oscillations were 
studied by Lindblom and Mendell jlCj by constructing an energy functional and 
computing the time-scales associated with the dissipative terms. The oscillation 
modes of superfluid analogs of the classical Maclaurin, Jacobi and Roche ellip- 
soids were derived recently by Sedrakian and Wasserman [fll| within the tensor 
virial method. The oscillation modes of superfluid ellipsoids separate into two 
classes corresponding to relative and co-moving (or center-of-mass, hereafter re- 
ferred as CM) motions of two fluids. The CM oscillations are identical to the 
oscillations of single-fluid ellipsoids and are undamped if one ignores the viscos- 
ity of the normal fluid. The mutual friction contributes only to the damping of 
the relative oscillation modes. One important feature of the latter modes is that 
they do not emit gravitational radiation as there is no mass transport associated 
with them. Our discussion of the tensor virial method is based on Ref. [ pd| . 

This review is organized as follows. In the remainder of the Introduction we 
motivate the averaged superfluid hydrodynamics and identify the relevant scales 
in the problem. In Sect. 2 we give a tutorial introduction to the Newtonian su- 
perfluid hydrodynamics, which is mainly built on the work of Bekarevich and 
Khalatnikov Sect. 3 introduces the tensor virial method and illustrates its 
applications to the superfluid ellipsoidal figures by computing the perturbations 
of the new terms in virial equations due to their two-fluid nature. We discuss in 
Sect. 4 the oscillations of superfluid Maclaurin spheroids including the effects of 
mutual friction and viscosity of the normal fluid. Sect. 5 is a brief summary. We 
refer the reader to the accompanying article by B. Carter for a review of relativis- 
tic models and an overview of the state of the art of the theory of superfluidity 
in neutron stars. 

1.1 Characteristic length scales 

The physics of neutron star superfluidity unfolds on a hierarchy of three distinct 
length-scales. The separation of these scales is useful, as often the physics of a 
neighboring scale enters a theory at a given scale in the form of phenomeno- 
logical constants, which can potentially be fixed by comparison with measured 
observables. 

At the microscopic level the physical scale of the order of fermi (fm= 10 - 
cm) is set by the nuclear forces. The long-range attractive interaction between 
nucleons leads to an instability of the normal state of the nuclear matter against 
Cooper pairing, in analogy to the microscopic Bardeen-Cooper-Schrieffer the- 
ory of superconductivity of metals. The nuclear forces control the size of the 
"elementary bosons" of the theory - the Cooper pairs, which appear as weakly 
bound states of two fermions near the Fermi surface. The size of a Cooper pair, 
the coherence length £, is of the order of 10 fm. It sets, obviously, the lower scale 
on which the hydrodynamic description of superfluids breaks down. On length 
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scales larger than £, the condensate of Cooper-pairs can be described by a single 
wave function "0( x )i i-e. , the condensate forms a macroscopically coherent state. 

At the local hydrodynamic level the relevant physical scale is set by the size 
of vortices - macroscopic quantum objects, whose fundamental property is the 
quantization of the circulation around a path encircling their core. The circu- 
lation is quantized in units of 2nh since the condensate wave-function must be 
single valued function at each point of the condensate. On writing ip = ipoe 1 *, 
the gauge invariant superfluid velocities can be expressed through the gradient 
of the phase of superfluid order parameter \ an d the vector potential, A 

V. = ^V Xr - (1) 

lm T m T c 

where e T = (e, 0) is the electric charge of protons and neutrons respectively, m T 
is their bare mass; r £ {n,p}, where n and p stand for neutrons and protons 
respectively. Applying the curl operator to Eq. ([!]) and implementing quantiza- 
tion of the circulation (the phase of the superfluid order parameter changes by 
2w around a closed path) we find 

curlv r = — v T V(5 (2) (x - x T? -) - -^—B = uj t , (2) 

3 

where Trh/m T is the circulation quantum, v T = lu t /uj t is a unit vector along the 
vortex lines, x TJ - defines the position of a vortex line in the plane orthogonal 
to the vector v T , 8^ is a two-dimensional delta function in this plane and 
B = curl A is the magnetic field induction. The j-summation is over the sites of 
vortex lines. Note that Eq. (||) treats the vortex cores as singularities in the plane 
orthogonal to v T , which is justified on scales larger than the coherence length of a 
condensate. The crossover from the local hydrodynamic scale to the microscopic 
scale can be studied within the Ginzburg-Landau theory as we briefly discuss 
below. 

For a single vortex the integral of Eq. @ completely determines the su- 
perfluid pattern; as this equation is linear, for larger number of vortices the 
superfluid pattern is a superposition of the flows induced by each vortex. The 
resulting net flow, obviously, depends on the arrangement of the vortices. It 
turns out that the integral of Eq. (||) on the local hydrodynamic scale is radi- 
cally different from the one on the scales involving larger number of vortices in 
a rotating superfluid. 

To appreciate the difference in superfluid patterns on different scales let us 
look for vortex solutions in a neutral condensate in cylindrical geometry. The 
condensate wave function has the form ip(x) = f{r)e ie in the polar-cylindrical 
coordinates (r,9,z); the neutron superfluid velocity upon integrating Eq. (^|) 
becomes 



(3) 
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The divergence of the superfluid velocity when r — > is avoided by introducing 
cut-off on scales of the order of the coherence length. The cut-off scale, for Fermi- 
liquids, can be understood by noting that an increase of v„ when r — > causes 
an increase of the kinetic energy of Cooper pairs which eventually becomes larger 
than the binding energy of a pair. The broken pairs will perform a rigid-body 
rotation with an angular velocity which scales as Cr and is regular in the r — > 
limit (here C is a constant). This crossover can be seen from the well-known 
solutions of the Ginzburg-Landau equation for the amplitude f(r): 



d 2 f 1 df_ _ 1_ 



where C, = r/£ n and is the size of the vortex core. The asymptotic solutions 
of Eq. (||) are 

/ ^=\l-(2C 2 )- 1 C-oo, (5) 

while numerical solutions for intermediate values of £ show that the conden- 
sate wave function is at half of its value in a homogeneous condensate when 
C = 1- Note the long-range nature of the superfluid vortex velocity, and the re- 
sulting slow fall-off of the density perturbation in the condensate. This behavior 
is specific to neutral condensates; for charged condensates the super-current is 
screened exponentially on length scales of the order of the penetration depth A. 
The solution of Eq. @ for a charged condensate is 



"p 



*i(tH ( 6 ) 



2m p X V A 



where K\ is the Bessel function of imaginary argument; as for r A, K\(r /X) ~ 
exp(— r/A), therefore the superfluid circulation decays exponentially; in the op- 
posite limit r <C A, K\{r/X) ~ X/r and Eq. (||) assumes a form identical to Eq. 

©■ 

At the global hydrodynamic level the relevant scales are of the order of the 
size of the system, which in neutron stars is of the order of kilometers. On 
these scales the hydrodynamic and thermodynamic variables are course-grained 
quantities, i.e. they are averages over a large number of vortices. The solution (||) 
does not minimize the energy E — L • f2 of a rotating superfluid, where E and L 
represent the kinetic energy and the angular momentum respectively. The energy 
acquires its minimum for a superfluid flow which to a high precision mimics a 
rigid body rotation i.e., v n — Qr, where r is the distance from the rotation 
axis; (small deviation occur only at the bounding surface of the superfluid). On 
the global hydrodynamic scales a transition to a continuum vortex distribution 
can be carried out on the right-hand side of Eq. (||) by defining vortex densities 
n T = (5^ 2 )(x — X-rj). Since for rigid-body rotations the curl of v„ is simply 
2/2, the number density of vortices in the neutron superfluid is related to the 
macroscopic angular velocity of the neutron condensate by the familiar Feynman 
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formula 

2m„i7 

nn = — fr . (7) 

For typical pulsar periods, 0.05 < P < 0.5 s, n„ ~ 6.3. x 10 3 P^ 1 — 10 4 -10 5 per 
cm 2 . The minute difference between the superfluid and normal angular velocities 
in a neutron star decelerating under external breaking torques is neglected here. 
For a charged superfluid Eq. (J2J) can be transformed to a contour integral over 
a path where v p = 0, as the super-current is screened beyond the magnetic field 
penetration depth A. Again, going over to the continuum vortex limit we find 

n p = — ~ 5 x 10 18 cm" 2 . (8) 

Note that the number of proton vortices per neutron vortex is n p /n n ~ 10 13 — 
10 14 independent of their arrangement. The energy of a bundle of neutron or 
proton vortices is minimized by a triangular lattice with a unit cell area 

n" 1 = (x/3/2) d 2 T . 

The length of a "basis vector" of such a lattice in a neutron condensate (the 
neutron inter-vortex distance) is 

nti 



For the inter-vortex distance in the proton condensate we find 

2i>„ N 1/2 



+ • (10) 

where B is the mean magnetic field induction. Using the estimates given in 
Eqs. (Q) and (|^) we find that the neutron and proton inter-vortex distances 
are d n ~ 10~ 2 — 10~ 3 cm and d p ~ 10 -9 cm respectively. For typical values 
of the microscopic parameters the penetration depth is of the order of 100 fm, 
therefore the conditions <C d n and £ p <C min(A, d p ) are satisfied and the use 
of the hydrodynamics on the local scale is valid. It is also clear that the global 
hydrodynamics can be applied on the scales that are much larger than d n (a 
fraction of millimeter). 

The remainder of this review concentrates on the physics of the global hy- 
drodynamic scale and on neutral superfluids only, as the dominant fraction of 
the moment of inertia of a neutron star resides in the neutron fluid and it plays 
the main role in the hydrodynamic oscillations of the star. Charged superfluids 
will be absorbed in the normal fluid of the star formally, as they are coupled 
to the electron liquid via electro-magnetic forces on short time-scales. Note that 
their role is crucial in controlling the mutual friction on the local hydrodynamic 
scale; however the physics of this scale will enter the theory on the global scale 
via phenomenological constants, which we will treat as free parameters of the 
theory. 
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2 Two-fluid Newtonian superfluid hydrodynamics 

The superfluid phases in neutron stars coexist with normal fluids whose inter- 
action with superfluid vortices leads to the effect of the mutual friction between 
a superfluid and normal fluid. A phenomenological description of this effect is 
based on the two- fluid dissipative hydrodynamics. A particularly simple and 
transparent formulation which, with some care, can be taken over to describe 
superfluids in neutron stars, was developed by Bekarevich and Khalatnikov for 
liquid He 4 ]T2| . It is interesting that the general form of mutual friction forces 
can be obtained by utilizing only the conservation laws and some reasonable as- 
sumptions on the form of the dissipation. Although the superfluids in a neutron 
star can be assumed to be at zero temperature, i.e., the number of quasi-particle 
excitations is small, they coexist with a normal liquid of electrons in the core 
and a nuclear lattice in the crusts. Hence, entropy is irreversibly produced due 
to various dissipative mechanism in the normal fluid even though the superfluid 
matter is effectively at zero temperature. Formally, the electron liquid and the 
nuclear lattice in the crusts take over the role of quasi-particle excitations in the 
superfluid hydrodynamics of liquid He 4 . 

The conservation of the combined mass of the two fluids is given by the 
continuity equation 

l+V.J-O, (11) 

where the net mass p = ps + Pn is the sum of the masses of superfluid and 
normal fluid, J — pn^n + Ps^s is the mass current (hereafter the indexes S 
and N refer to superfluid and normal fluid, respectively.) The total momentum 
conservation is 

where Pik is the stress energy tensor. The time evolution of the entropy, S, of 
normal fluid can be written as 

dS _ _ R . , 

- + V.*v„ = -, (13) 

where R is the dissipative function, T is the temperature; finally, the conservation 
of the energy, E, reads 

^ + V-Q = 0, (14) 

where Q is the energy current. Equations above should be supplemented by the 
Euler equations for the superfluid and the normal fluid 



Ps 

PN 



(v s ■ V) • v s 



dt 
dt 



+ (vat • V) • Vat 



-^Vp-p s V^ + F, (15) 
P 

■— Vp - pjvV^ + jyjv^vjv - F, (16) 
P 
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where F is the mutual friction force, and rj^ is the viscosity of the normal fluid, 
4> is the Newtonian gravitational potential. To determine the unknowns in the 
hydrodynamic equations, let us write the total energy of the fluid in the frame 
in which the normal fluid is at rest as 

E=^pv 2 s + (J-pv s )-vs+£, (17) 

where the internal energy £ is given by the second law of thermodynamics as 

d£ = TdS + pdp + (vjv - v s ) • d(J - pv s ) + Adw. (18) 

The energy due to the vorticity is represented by the term which is proportional 
to ui = V x vs- Differentiating Eq. ( |l7|) with respect to time and eliminating the 
time derivatives using the conservation laws above we recover the conservation 
of the energy 

+ ( J - pv N + V x Av) ■ {F + [(v s - v w ) x w]} = 0, (19) 

where P' ik is the part of the stress tensor associated with the vorticity; the 
explicit form of the energy current is not indicated since it will not be used in the 
following. Since the dissipative function R must be positive, the remaining terms 
on the right-hand side of Eq. ( {l9| ) must be quadratic forms for small deviations 
from equilibrium. This implies that the most general form of the mutual friction 
force is 

F = - [uj x (V x Av)] — P[v x [u x (vat - v s - V x Av) ] ] 

- ft [w x ( VA r - v 5 - V x Av)] + (}" v ■[w(v JV -v s -Vx Av)] , (20) 

where (3, (3' and [3" arc phcnomcnological coefficients. On substituting the mutual 
friction force in the Euler equation for the superfluid, Eq. (|l5|), we see that the 
vorticity propagates with a velocity v^, that is 

|^ = Vx(v L x^), (21) 

which is defined, assuming f3" <C (3 ,f3', as 

Vi = v s + V x Av + (3'(w N — v s — V x Av) 

+ (3[u x {v N - v s - V x Av)] . (22) 

The latter equation can be put in a form reflecting the balance of forces acting 
on a vortex 

Ps [(vs + V x Av - v L ) x lu] - 77 (v L - vat) + 7/ [(v L - vat) Xu]=0, (23) 
with the new phenomenological coefficients rj and rj defined as 

p = vps^ &l = 1 _ P8»{pa»-4) (24) 

rj 1 + (p S LO - T)') rj 1 + (p S LO - r)') 
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The first term in Eq. (23) is a non-dissipative lifting force due to a superfiow past 
the vortex (the Magnus force). The remaining terms reflect the friction between 
the vortex and the normal fluid. The coefficients rj and rj', therefore, measure 
the friction parallel and orthogonal to the vortex motion in the plane orthogonal 
to the average direction of the vorticity. A nonzero (3" implies friction along the 
average direction of the vorticity, which is possible if vortices are oscillating, or 
are subject to other deformations in the plane orthogonal to the rotation. One 
may assume, at least under stationary conditions, that f3" <C /?,/?'. 

In the both limits of either strong coupling (77 3> psk>) or weak coupling 
(77 <C ps^ 1 ) between a vortex and the normal fluid, one finds that (3 — > as a 
function of rj, with the maximum /3 max = 0.5 at rj = ps^. In the strong coupling 
limit (3'irj) — ► 1, while in the opposite weak coupling limit (3'(rj) — > ( generally 
we assume that the quasi-particle-vortex scattering kinematics implies rj' <C rj 
and that for the relevant densities rj' <C pslo). Note that f3'{rj) approaches its 
asymptotic strong-coupling values quadratically, while (3(jj) does so linearly; the 
asymptotic behavior for large 77's, therefore, is dominated by /3(rj). 



3 Virial equations and their perturbations 

Virial equations of various order are constructed by taking moments of the hy- 
drodynamic equations. Since the computation of the perturbations of the these 
virial equations is central to the theory of superfluid ellipsoids we review this 
somewhat technical issue in this section. The reader who is interested only in 
the physics of superfluid oscillations can proceed to the next section where we 
discuss the oscillations of Maclaurin spheroids. 

The equations of motion ( |l5| ) and (|l6|), written in a frame rotating with 
angular velocity u> relative to some inertial coordinate reference system, can be 
combined in a single equation 



d d \ dp a dcf) 1 8\ujxx 



2 



dt ' dxj J dxi ^ a dxi 2^ a dxi 

-(- ^Pa^ilm u a ,iO m + Fa^i, (25) 

where the subscript a € {S, N} identifies the fluid component, and Latin sub- 
scripts denote coordinate directions; p ai p a , and u Q are the density, pressure, and 
velocity of fluid a. The two fluids are coupled to one another by mutual gravita- 
tional attraction and the mutual friction force F Q ^ [Eq. (20)]. The gravitational 
potential 4> is derived from 

V 2 = V 2 (0 S + M = 4^G[p s (x) + Mx)]; (26) 

the individual fluid potentials 4> a obey V 2 </> Q = 4irGp a . For a normal-superfluid 
mixture the mutual friction force, written in components, is 



FsN,i = -Psusf3ij{us,3 - u N j), 



(27) 
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where the mutual friction tensor is 

Pij = PStj + (3'e ijm u m + (/?" - flvivj, (28) 

with /3, j3' and f3" being the mutual friction coefficients, and u>s = vcus = 
Vxu s . 

The Euler equation ( p5| ) can be extended to include external gravitational 
sources, for example the tidal potential of an external point source of gravity 
acting on an ellipsoid (Roche ellipsoids). We will not discuss here the stability 
and oscillations of superfluid counterparts of the classical Roche binaries; their 
relative oscillation modes are derived in Ref. jo]. 

The general strategy for finding the equilibrium shapes of ellipsoidal figures 
and modes of their oscillations within the tensor virial method consist of (i) 
constructing moments of the hydrodynamic equations describing fluid motions 
in the rotating frame; (ii) computing Eulerian perturbations of the resulting 
virial equations; (iii) expressing these perturbations in terms of the virials of 
various order; these are defined as the moments of the Lagrangian displacement 
£ ai of fluid a: 

V a ,i = / d*xp£ u (first order) (29) 
Jv a 

V a ,i;j = / d 3 xp£ ai Xj, (second order) (30) 

V aA .j k = / d 3 x pi al x 3 x k , (third order) (31) 
Jv a 

The advantage of using the homogeneous ellipsoidal approximation is that the 
perturbations of the gravitational energy tensor of an ellipsoid can be expressed 
in terms of the index symbols defined as (cf. EFE Chap. 3) 

Aijk... — aia 2 a 3 / — 2 — 2 — — 2 — , (32) 

7o A(af + u)(aj + u){a 2 k + u) . . . 

f°° udu . . 

tiijk... = a!a 2 a 3 / — — — — , (33) 

Jo A(aj + u)(aj + u){a 2 k + u) . . . 

where A 2 — (a\ + u){a\ + u){a\ + u) and oi, 02, and 03 are the semi-axis of 
the ellipsoid. This strategy is described in detail in EFE. Below, we concentrate 
on its extension to superfluids with an emphasis on the new effects of mutual 
friction and mutual gravitational attraction of the superfluid and normal fluid. 

3.1 First order virial equations 

On taking the zeroth moment of Eq. (|25| ) which amounts to integrating over V a 
we obtain the "first order 'virial' equation" 2 

(X„ d3x PaUa '*) = 2£lim ^ m J d 3 xp a u a .i + (Q 2 5 lj - QiQj) J d 3 xp a xj 
2 The word "virial" is in quotes because the equations are intrinsically dissipative. 
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-(1-tfa/j) / d 3 x Pa ^+ f d 3 xF a p. l . (34) 

Note that we impose the boundary condition p a = for each fluid on the bound- 
ing surface of V a . The fluids are not restricted to occupy the same volume. Apart 
from simple doubling of the number of the incrtial forces, which do not couple 
the two fluids, there are two forces that do couple them: gravity and friction. The 
net mutual gravitational force between the fluids vanishes only if they (i) occupy 
the same volume and (ii) have densities that are proportional to one another (i.e. 
ps oc pn). The mutual friction force is nonzero as long as the fluids move relative 
to one another. Although the mutual friction force is nonzero only in the overlap 
volume of the two fluids - a restriction which is necessary to derive conservation 
of total momentum for the combined fluids - it would be effective throughout 
the entire volume of fluids because the force is mediated by a macroscopically 
extended vortex lattice. 

For isolated single-fluid ellipsoids the first harmonic oscillations are trivial, 
since they correspond to the motions of the center-of-mass of an ellipsoid and 
can be eliminated by a transformation to the reference frame whose origin is 
located at the center-of-mass of the ellipsoid. For two-fluid ellipsoids the fluid 
motions include the counter-phase (relative) oscillations of two fluids, which can 
not be eliminated by any transformation. These are the only new type of oscil- 
lations for the superfluid Maclaurin and Jacobi ellipsoids (the solitary ellipsoids 
with vanishing internal motions). In the case of ellipsoids in an external grav- 
itational field (e.g. the Roche ellipsoids), the CM motions are not trivial any 
more, since the external (inhomogeneous by assumption) source of gravitational 
field breaks the translational symmetry of the problem. Hence, apart from the 
relative oscillations of two fluids, the CM oscillations become non-trivial. 

Consider first the variation of the first order virial equation under the influ- 
ence of perturbations. The variations of the inertial terms proceeds in full analogy 
to EFE. Here we concentrate on variations of the new terms corresponding to 
the mutual gravitational attraction and mutual friction. The variation of the 
first force is 

d(t>f 3 f 3 . . . . d f d 3 x'pp(x.')(xi - Xi) 



i f d 3 xp a ^ = - [ d 3 xp a (x.)Ui(^)^- I 
Jv a oxi J Va dxi J v 



f * / \* r \ ® f d 3 x'p a (x')(xi - x'A ,„„, 

which is manifestly antisymmetric on a <-> (3. Assuming V a — Vp — V and 
Pa = /qP( x ) m the background equilibrium, we can simplify this to 



6 l v d3xPa ^ = fafp //Mx)KwW-( tt ,i(x) 



V 

d_ f d 2 x'p(x , )(x i ~x' l ) 
i Jv 



dxi J v Ix-x'l 3 



(36) 



Although we simplified the final answer by assuming that the fluids occupy 
identical volumes and have proportional densities in the background state, we 
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could not have derived the correct perturbation of the first order virial theorem 
if we had not allowed the volumes to differ. 

For uniform ellipsoids, we can simplify the mutual gravitational term further. 
First, we note that the second integral is simply the derivative of the gravitational 
potential which at any interior point of a homogeneous ellipsoid is 



k=l 



x) = -irGp(l-^A k 4), (37) 



where I is a constant. Thus, the mutual gravitational contribution to the equa- 
tion of motion for the perturbed center-of-mass is 

2nGp 2 A l f a f p [ d 3 x(Z 0ti -Z a ,i). (38) 
Jv 

If in the background state, the two fluids move together or are stationary, the 
variation of the mutual friction force becomes 



5 J d i xF a i3. i = -S a p^ f s J d 3 xp(x)Ljsf3ij(£sj - £nj) 



(39) 



where S a p = if a = /3, 1 if a = S and (3 — N, and —1 if a — N and (3 = S. 
Collecting terms we find the first order virial equation 



fa — 2(zH m f a f2 m . Vq:,Z "1" (^ &ij ^i^2j^)faVa,j 

- 27rGpAJ a f (3 (V aii - V^) - S a pf s us(3ii {V a .j - V PJ ) . (40) 

The CM and relative motions can be decoupled by defining 

Vi = f s V s ,i + f N V N , h Ui = V s ,i-V N ,i. (41) 

The CM motions of two fluids are trivial (as they can be eliminated by a linear 
transformation of the reference frame) and, therefore, Vi = 0. The virial equation 
describing the relative motions of the two fluids is 

^Ui = 2e llm n m jUi + (I2 2 <% - Q i Q j )U j - 2A i XJ i - 2Q^ i ~U j . (42) 

From the latter equation it is straightforward to compute the first harmonic 
relative oscillation modes of irrotational ellipsoids, which is done in the next 
section for Maclaurin spheroids. 

3.2 Second order virial equations 

Taking the first moment of Eq. ( p5|) results in the second order 'virial' equation 

+ ^T a ,ij + 8ijH a + Ma.ij + (1 — 5 a p)M a p,ij + Fafidj, 

(43) 
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where 



Ia,ij — / d X p a X^Xj (44) 

n a = f d 3 x Pa (45) 
Jv a 
I r 

Ta.,ij = 7) / ^ % PaV , a,i'Ua.j (46) 
Z JV a 

G f d 3 xd 3 x'p a (x)p a (x')( Xi -x'Jixj - x'j) 
■M a ,ij =~i7 1 .,3 ( 47 ) 



M 



d 3 x'p a (x)pp(x / )x J (x l - a;-) 

v a Jv (i 

Fap.ij = I d 3 xXjF aj 3 A . (49) 



_ r f * f d 3 x'p a (x)p fj (x / )x :j (x l -x' i ) 

L L ix-x'|3 (48) 



af3,ij 



When there is only one fluid present, this equation reduces to the one found 
in Chap. 2 of EFE. Again we consider only variations of the new terms in the 
second order virial equation due to the mutual gravitational attraction (A4 a p,ij) 
and mutual friction {T a p,ij)- The first variation is 

d f d 3 x'p(-x.')(x i -x' i )(x j -x' j ) 



SM af 3,ij = -Gf a fp^J^d 3 xp(x)£ a (x)-^- J 



I ./i |x-x'| 3 



+ y/ 3 ^(x)[C(x)-^(x)]^ J i),(50) 

where we have specialized to backgrounds with proportional densities and iden- 
tical bounding volumes. The first term in the brackets can be combined with 
SM a,ij- The resulting equation can be written more compactly in terms of the 
functions 

d 3 x'p(x)(x l - x'^Xj - x'j) 



dxi 

which are related by 



b, ,j y x x /,„ (5i) 

8Dj G r d 3 x'p{x')x' J {x l -x' i ) 
Jv |x-x'| 3 



?El- b ..- x .°* 
d Xl ~ 13 3 dx. 



$ij %j ! (53) 



wc find 



'•Ma.i 3 + (1 - 5 a p)5M a p.ij = -f a J d 3 xp£ a j^^- 



+ Up [ d 3 xp(Ui ~ fai)S^r-- (54) 
Jv dxidxi 
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For the uniform ellipsoids (cf. EFE, Chap. 3, Eqs. [125] and [126]), 



irGp 
By 

TlGp 

■jrGp dxidxi 



^ fe=i ' 
2BijXiXj + <5y ( A, — ^ ' A^x^ j . 

2B i] {5ux ] + SjiXi) - 2a i 5ijA i ixi. 



(55) 

(56) 
(57) 



Using these results, and defining symmetric in their indexes second order virials 

as 



we finally obtain 

SMa.ij + (1 - S a p)6M al 3,ij 



(58) 



7TG/9 



fa ( 2BijV a .ij — a { Sij 



1=1 



-Ojfafp 



2Aij(V a ,ij - Vj3,ij) + Si 3 ^2 Aii{V a ,u - Vp.ii) 



For the perturbations of mutual friction force we find 



1 



S I d 3 xxjF af3ti 



Sapfs / d 3 xpuj s xj/3tk 
Jv 



dt dt 



(59) 



(60) 



For perturbations of uniform ellipsoids, u>s and p are independent of position in 
the unperturbed background, and we may also assume that /3ij is constant; for 
backgrounds in which there are no fluid motions 



S / d 6 XXjF al3 ^ = -SapfsP^sPik 

JV a 



dt 



dt 



(61) 



Thus the second order virial equation for a fluid a, in a frame rotating with an 
angular velocity f2, is 



dt 2 



dt 



^ i=i ' 



2A ij (V a , ij - V^ij) + 5 tl Aii{V a ,u - Vfi,ii) 



- OjfafpwGp 

- S a fjf a UJsPik-^[ la.t:j - Vl.l.-.J )■ 



(62) 
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We can replace these equations with a different set by defining 

Vi-j = fsVs,i ; j + /nVn^-j Ui-j = Vs,i-j — Vpf t i-j. (63) 
In terms of these new quantities we find 

^2 = ^timftm ^' + ft Vij — QiQkVkj + SijSII 

- nGp(2B ij V ij - a% 3 £ ^VjA , (64) 



i=i 



N 



<i-i ,:j .. n du v , 2 (&n s sn 

— ,, 9 — &£ilm**m—r. r J/ t/jj — llillll k j + Ojj — j j — 

«e at V Js JN 

- 2nGpA i U ij - 2Qp ik jU k , r (65) 

The first equation is identical to the second order virial equation for a normal 
inviscid fluid. The second equation is specific to superfluids and contains all the 
new modes of relative oscillations between the normal fluid and superfluid. It is 
clear that the separation of the oscillation modes in the CM and relative modes is 
the result of the symmetry of the two- fluid hydrodynamic equations with respect 
to the interchange a <-> /3. If this symmetry is broken the two classes of modes 
mix. We shall consider below the effect of the viscosity of normal fluid which 
by definition acts only in the normal component and hence breaks the a <-> (3 
symmetry. 

The second order virial equation for viscous fluids, quite generally, acquires 
the term 

-n f v a v - x ( dUm , dUa - k 2 dUa > 1 x \ IRR\ 
r a ,tj = J Pa.ijdxi, P a . ik = b a , N p N v y-Q^ + —q^ 3~dxi~ ) ' 

which is called the shear-energy tensor; v is the kinematic viscosity^. For back- 
ground states which are stationary and without internal motions the variation 
of the stress-energy tensor is 

SV^ = 6 a , N J^ PNV - ^_ + _ - . (67) 

It is impossible in general to express the variations of the stress-energy tensor in 
terms of the virials V a i-j. However, in the low Reynolds- number approximation, 
this tensor can be approximated in a perturbative manner using as the leading 
order approximation the proper solutions for the displacements corresponding 
to the inviscid limit. Since the latter are linear functions of the virials, £jvi = 



3 We use the same symbol v for the kinematic viscosity and for the unit vector along 
the vortex circulation; no confusion should arises as the latter quantity does not 
appear in the virial equations. 
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X)fc=i ^N,i-,kXk/M]^a\, with Mjv being the mass in the normal fluid, one finds 
8V atij = 5uS a>N ± fei + ^f-^j (68) 

in the incompressible limit. Thus, the second order virial equation which includes 
the viscosity of the normal matter becomes 

fa = ^il 7n ^ 7n f a ^ ~ ^ faVa.ij ^i^kfaVa.kj ~t~ 3ij3n a 



1=1 



a? a? 



(69) 



Apart from the last term, the remaining terms in Eq. ( p9|) manifestly preserve 
the symmetry with respect to the interchange a <-> (3; note that they might have 
different parities under this transformation. The last term breaks this symmetry 
as the viscosity acts only in the normal fluid. 

3.3 Third order virial equations 

To obtain the third order virial equation we take the second moment of Eq. ( p5| ) 
and integrate over V a : 



dt 



where 



d Xp a U a iXjXk J — 2^-ilra^m ^ -^Pa^a,lXjXk j ~\~ ^ la.ijk 

— f2if2il a jjk + 2(T a ,ij-k + T a ,ik;j) + Sijlla^k + SikII a j 
+ M-aP.ijk + 3~af3,ijk, (70) 



Ia,ijk — / d X p a XiXjX]c (71) 

Jv a 



n a ,i = / fspa.i (72) 

JV a 

I r 

1~a,ij;k — 77 / ^ ^Pa^a.i^a.j-^k (73) 



l'„ 



2 

M a pjjk = - I d^XXjXkPa^r 1 - (74) 
^aff,ijk = / d 3 X Fap.iXjXk- (75) 

Jv a 
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Below, we compute only the perturbations of the tensors in the last line of Eq. 
(fro|), which correspond to the gravitational potential energy and the mutual 
friction; the perturbations of the remainder terms is computed in analogy to 
EFE. 

For the Eulerian perturbation of the gravitational potential tensor we have 



5 I d xXjXkPo 



dxi 



SG / d 3 XXjXkPa(x) 



'13 



+ / d 3 x'pp{x')^ <] 



v 



'13 



(76) 



Assuming V a = V and p a = f a p(x) in the background equilibrium, and defining 
[cf. EFE, Chap. 2, Eqs. (14) and (22)] 



Bij = G / cfVp(x') 



,/|3 



V lk .j = G [ d?x'p(x')x 
Jv 



X — X' 

, (Xi — x'^lxk — x' k ) 



r l\3 



(77) 
(78) 



one finds for the 'self-interaction' term that 



-5 



XX j Xfcp G 



a lty = J d3xp ^ a ' l lix l ( Bi i Xk+T>i r< k ") 

- l -H d 3 xp(x)£ a j^- (BikXj + V ik .j) = fl(5M a ^ k + 5M a , ik]j ). (79) 

To arrive at the symmetric in the indexes k,j we used the identity 

BijX k + V lk .j = BikXj + T>ij. k . (80) 

The perturbation of mutual interaction term in Eq. ( f76| ) , assuming again V a = 
Vp = V and p a — f a p(x) in the background equilibrium is 



,3 

a xxjXfcp a 

OXi 



fafp I d 3 xp(x)£ al (x) — {BijX k + T>ik-,j) 



dxi 



fafp I d 3 X P (x) [C a ,,(x) " ^(X)] 

Jv 



d 2 v 



jk 



dxidxi 



where 



V 



jk 



G I 



(81) 



(82) 



Combining Eqs. © and © we find 

-8 j d 3 XXjX k Pa^- = - fa ( d S Xp(x)£ a (BtjXk + T> ik -j) 

JVr, " x i JV ' " x l 



Up f d 3 x P {x) [^(x)-^(x)] 
Jv 



dxidxi 



(83) 
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For the perturbation of the mutual friction tensor, assuming stationary back- 
ground equilibrium, we find 



S^aBAjk = S I d 3 xF a B ,XjX k 



Safifs / d XXjX k p(yL)L0sf3u 



NJ 



dt dt 

Putting together all the terms we arrive at the third order virial equation 



(84) 



J°> j79 = ^ilmiimja h Oijdll a ^ k + Oi k OU a j 



dt 2 



dt 



f d 

+ {Q 2 5 a - f2iOi)f a V a ,ijk ~ fa J d 3 xp(x)£ al — {Bi 3 x k + T>ik-,j) 

-Up I d 3 xp(x) [^(x)-^(x)] 
Jv 



dxidx. 



dVs,k]k dVN,l-jk 



(85) 



dt dt 

where the symmetric in its indexes third order virial is defined as 

*a,ijk — Va,i;jk ~t~ *ct : j\ki ~t~ ^a,fc;ij- 

To separate the CM and relative motions of the two fluids introduce the virials 

Vi-jk = fsVs.i;jk + fNVN,i-jk Ui-jk = Vs,i;jk — Vjy,i;jk- (87) 

The new set of equations in terms of these virials is 



d 2 Vj-jk 
dt 2 



dV, 



i-jk 



dt 



S tj Sn k + SikSIIi + (f2 2 S tl - QiQi)Vi jk 



f 9 

J d 3 xp(x.) [fs£ s ,l + InZnj] (BijX k + V ik .j) . 



(88) 



and 



d 2 Ui-jk 
dt 2 



2e am Q m — I 1 + -7— ) wsPu 

JN 



dUi-jk , r / Slls.k 5II N . k 



TN 



+ S, 



sn 



S,j 



sn 



{Q 2 5 a - f2 t t2i)Ui jk - 2-KGpAiU ijk 



fs In 

f d 
2 J d 3 xp( X ){£ s . l - £ N l ]—(B l3 x k + V ik;j ) . (89) 



To obtain the last term in Eq. ( p9[) we used the relations [cf. EFE, Chap. 2, Eqs. 
(29) and (28)] 



dT>j k dV k 



dxi 



dxi 



dV k dcf> 
— — = B ik + x k ——, 

axi oxi 



(90) 
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and the explicit expression for the gravitational potential of an ellipsoid, Eq. 
(|37|). The terms in the last lines of Eqs. ( p8| ) and ( |S9| ) can be worked out to a 
form involving a linear combination of virials and index symbols, however the 
present form already makes clear that they will involve the virials describing the 
CM and relative motions, respectively. 

4 Small amplitude oscillations of superfluid Maclaurin 
spheroid 

In this section, we specialize our discussion to Maclaurin spheroids, the equilib- 
rium figures of a self-gravitating fluid with two equal semi-major axis, say a\ and 
0,2, rotating uniformly about the third semi-major axis 03 (i.e. the X3 axis). For 
these figures in many cases analytical results are available. The superfluid oscil- 
lations of more complicated non-axisymmetric figures like the Jacobi and Roche 
ellipsoids require numerical analysis which is beyond the scope of this review 
(seeRef. Q). The sequence of quasi-equilibrium figures of Maclaurin spheroids 
can be parameterized by the eccentricity e 2 = 1 — a§/a 2 , with (squared) angular 
velocity fi 2 = 2e 2 B 13 , in units of (tt/oG) 1 / 2 . 

Surface deformations related to various modes can be classified by corre- 
sponding terms of the expansion in surface harmonics labeled by indexes l,m. 
We shall concentrate below on the first and second harmonic surface deforma- 
tions correspond to I = 1, 2 and — 1 < m < 1, — 2 < m<2 respectively. 

4.1 First order 

If the time-dependence of the Lagrangian displacements is of the form 

= tafa)**, (91) 

then the characteristic equation for the first order relative oscillation modes 
becomes 

X 2 U t = 2e llm Q m \Ui + (fPSij - Qif2j)Uj - 2A l U l - 2!2j3 ij \U j , (92) 

where all frequencies are measured in units (irGp) 1 ^ 2 , = ( 1 + /s/ / jv)/% j and, 
since we assumed no internal motions in the background equilibrium, = 2fl. 
The CM oscillations are trivial as they can be always eliminated by a transfor- 
mation to the center-of-mass reference frame. 

Assume that the ellipsoid is rotating about the X3 axis. Then, writing Eq. 
( p2| ) in components we find 

(A 2 + 2f2(3X + 2Ax - r? 2 ) Ui - 2f2(l - (3')XU 2 = 0, (93) 
(A 2 + 2f2(3X + 2A 2 - r? 2 ) U 2 + 212(1 - j3')XUi = 0, (94) 
(A 2 + 2Qp"X + 2A 3 ) U 3 = 0. (95) 
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The equations which are even and odd with respect to the index 3 decouple. For 
the perturbations along X3 the relative displacement vanishes, U3 = 0. Eq. ( p5| ) 
(which is odd in index 3) gives, on writing A = ia, 



CT? d 2 d = ±\J 2A 3 - P" 2 f2 2 + 0"H. (96) 

The first order odd parity oscillations are hence stable and damped if j3" < 
2 A3 /SI 2 . For Maclaurin spheroids 2 A3 /SI 2 > 5.040; the lower bound corresponds 
to eccentricity of the ellipsoid e = 0.865. For Jacobi ellipsoids this minimal value 
is slightly lower, 2A3/SI 2 = 4.148, and occurs at the point of the bifurcation of 
the Jacobi sequence from the Maclaurin sequence where the axis-ratio is defined 
by Cos _1 (a3/ai) = 54.48. Since the /^"-coefficient is the measure of friction along 
the average direction of the vorticity, it is reasonable to assume that $" <§; $, ft'; 
and since (3 < 1/2 and f3' < 1, we may conclude that the odd modes are stable, 
unless unphysical large values of j3" are adopted. The characteristic equations 
for the modes even in index 3 is 



A 4 + 4X2,3 A 3 + 2 



A 



(Ai + A 2 ) + Sl 2 {\ + 2/3 2 - 4/3' + 2/3' 2 ) 
+4/3/2 (At + A 2 - Q 2 ) A + {2A 1 - f2 2 )(2A 2 - Q 2 ) = 0. (97) 
For Maclaurin spheroids (A\ = A 2 ), upon writing A = ia, the solution becomes 



<2 cn = ipn± \J2A l - Q 2 [l - (3 2 - (1 - p') 2 }. (98) 

These modes represent stable, damped oscillations since the inequality 1 — /3 2 — 
(1 — /3') 2 < 2A\f SI 2 is always fulfilled. Indeed, the left-hand side is always larger 
than unity, while the maximal value of the right-hand side is 1/2. The latter 
upper limit is easy to deduce by minimizing the left-hand side of the inequality 
with respect to < 77/ ps^s < 00 defined via the relations (p4|). 

The analysis above shows the principal distinction between the modes gov- 
erned by the first and second order virial equations. At the first order the dis- 
sipation determines the unstable configuration (although it appears that this is 
never the case for realistic values of the friction coefficients). At the second order, 
as we shall see below, the dissipation only determines the rate of the instability, 
but not the stable configuration itself. 



4.2 Second order 

Second order harmonic deformations correspond to I = 2 with five distinct values 
of m, —2 < m < 2. Again, let us assume time-dependent Lagrangian displace- 
ments to have the form ([Jl]). The characteristic equation for the second order 
oscillation modes become 

^ 2 ^;i = 2eiim^mA + f2 2 Vij — Qif^kVkj 
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V fs f. 



N 




+ 5 tj SII - nGpUBijVij - aj8 VJ ^ A U V U 
^ 1=1 

_ .- )( , /vA ^ + ^ + hvf N f s \ ( ^J- + % ] , (99) 

. ( 5iis sn N 

+ o, 



2f20 ik XU k;j + 5i/A ^ + - 5^/sA + , (100) 



where the frequencies are measured in the units (71-/9G) 1 / 2 . Eqs. (99) and ( |100| ), 
which (if written in components) constitute a coupled set of 18 equations each, 
contain all the second harmonic modes of isolated, incompressible, and irrota- 
tional superfluid ellipsoids. In the next sections, we concentrate on solutions of 
these equations for the special case of Maclaurin spheroids, i.e. the case where 
the axial symmetry about the axis of rotation is assumed. 



4.3 Transverse shear modes (I = 2, m = |1|) 

These modes correspond to surface deformations with \m\ = 1 and represent 
relative shearing of the northern and southern hemispheres of the ellipsoid. They 



are determined by the eight components of the Eqs. (|99j) and ( 100 ) which are 
odd in index 3; i.e. V 3 -i, Vi- 3 , U 3 -i and Ui- 3 , where i = 1,2. The odd equations 
for the virials describing the CM-motions are 

(A 2 + f N vX + 2B 13 ) V 13 - (A 2 + -yf N v\ + 2B 13 ) V U3 

-f N f S vXU 13 + f N fsl^U h3 = 0, (101) 
(A 2 + f N vX + 2B 23 ) V 23 - (A 2 + 1 f N v\ + 2B 13 ) V 2 -, 3 

-fNfsv\U 23 + f N fsiyXU 2 , 3 = 0, (102) 
(A 2 - 7/a^A) V U3 - 2f2XV 2 . 3 + (2B 13 - Q 2 + f N v\) V 13 

+fNfsl^U 1 . 3 - f N fsvXU 13 = 0, (103) 
(A 2 - 1 f N v\) V 2 , 3 + 2nXV h3 + (2B 13 -Q 2 + f N vX) V 23 

+fNfsl^U 2 . 3 - f N fsPXU 23 = 0, (104) 

where 7 = l — a 2 /a 3 and we have redefined the kinematic viscosity as v' = hv/a\ 
and dropped the prime in the equations above. Relations ((5^) were used to 
manipulate the original equations to the form above. For the virials describing 
the relative motions the odd parity equations are 

(A 2 + 2QP"X + f s vX + 2A 3 ) C/13 - (A 2 + 2f2p"X + 7/s^A) C/ 1;3 

-vXV 13 + jy\Vi ]a = 0, (105) 
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(A 2 + 2/2/3" A + f s v\ + 2A 3 ) U 23 - (A 2 + 2(2(3" \ + 7/gi/A) U 2 , 3 

-^723+7^3,3 = 0, (106) 

(A 2 + 2Q~pX - yfsuXj U 1;a + {2A t ~ tt 2 + f s vX)U l3 - 212(1 - P')\U 2 ., 3 

+7i/AVi;3 - z^A7i3 = 0, (107) 

(A 2 + 2Q~pX - 7/s^A) U 2 ; 3 + (2Ai - f2 2 + f s vX)U m + 212(1 - P')XU 1;3 

+jv\V 2 ; 3 - vXV 23 = 0. (108) 
According to the symmetries of the original virial equation (|69|), the two sets 



(|01j)-(|0J) and fllOSD-fllOSD decouple in the limit v -> 0, as they should. The 
dissipation in the first set is driven by the viscosity of the normal matter; the 
superfluid contributes to the damping of the CM modes indirectly, via their 
coupling to the relative oscillation modes. In the second set the normal matter 
viscosity directly renormalizes the mutual friction damping time scale (212/3 — > 
2/2/3 — 7 /sf), thus reducing the damping of the relative modes. Note that this 
renormalization vanishes for a sphere since then 7 = 0. One important feature of 
each set, which remains preserved when they are coupled, is the balance between 
the tensors describing the perturbations of the rotational kinetic energy and 
gravitational energy; in the first set only the two- index symbols enter (-By), 
while in the second one appear only the one-index symbols (A,). As a results the 
neutral points (if any) along a sequence of ellipsoids (parametrized in terms of 
eccentricity) remain unaffected by the coupling between the different sets. This 
implies that as long as there are no neutral points for the relative transverse- 
shear modes in the uncoupled case, the conclusion about their stability can not 
be affected by the viscosity of the normal component. The CM modes do not 
show neutral points along the Maclaurin sequence and therefore their stability 
is guaranteed. 

In the absence of the viscosity the components of Eq. (|99|), which are odd 
in index 3, decouple into two separate sets. The first set for virials Vij, which 
describes the CM motions of the fluids is identical to the one found in EFE: 

(A 2 + 2£?i 3 ) V13 — (A 2 + 2Bi 3 ) Vi-3 -0, (109) 
(X 2 + 2B 23 )V 23 - (X 2 + 2B 13 )V 2 , 3 =0, (110) 

x 2 v U3 - 2nxv 2 , 3 + (2B 13 - n 2 ) v 13 = o, (in) 

X 2 V 2 , 3 + 20XV h 3 + (2B 13 - n 2 ) V 23 = 0. (112) 

The corresponding modes are described in EFE (see also Fig. 1 below). The 
second set, which describes the relative oscillations of the fluids, is 

(A 2 + 2Qj3"X + 2A 3 ) Uia - (a 2 + 212/3" a) U 1;3 =0, (113) 

(A 2 + 2Qj3"X + 2A 3 ) U 23 - (A 2 + 212/3" A) U 2 , 3 =0, (114) 

(A 2 + 2Q~px) U x .3 + (2A 1 - f2 2 )U 13 - 212(1 - P')XU 2 , 3 = 0, (115) 



The tensor virial method 23 




Fig. 1. The real (upper panel) and imaginary (lower panel) parts of the CM and 
relative transverse-shear modes of superfluid Maclaurin spheroid as a function of ec- 
centricity for values of rj/cosps = 0.0 (solid line), 0.5 (long-dashed line), 1 (dashed line), 
and 50 (dashed- dotted line). The left panel corresponds to the solutions in the inviscid 
limit v = 0. The right panel corresponds to the case where v — 4/3f2. The fraction of 
the normal fluid /jv = 0.1. 

(A 2 + 2i7/3A) U 2 -s + (2Ai - n 2 )U 23 + 2(2(1 - p')\U 1;3 = 0. (116) 

Let us concentrate first on the second set and consider the limit of zero mu- 
tual friction (i.e. the case where the two fluids are coupled only by their mutual 
gravitational attraction). The characteristic equation can be factorized by sub- 
stituting A = icr to find 

a [a 2 - 2(A 1 + A 3 ) + f} 2 } ± 2f2(<7 2 - 2A 3 ) = 0. (117) 

The purely rotational mode a = fl decouples only in the spherical symmetric 
limit where A\ = A 3 . If only axial symmetry is imposed then the third order 
characteristic equation is 



a 3 ± 2n<j 2 + [-2 (A 1 + A 3 ) + n 2 ] a T 4 A 3 fl = 0. (118) 
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It is easy to prove that the modes are always real. Three complementary modes 
follow from Eq. ( |118| ) via the replacement J? — * — Q. Fig. 1 shows the real and 
imaginary parts of the transverse-shear modes along the Maclaurin sequence. 
The left panel shows the results when v = 0. In that case the relative modes, 
which start for fl — > at 1.63, are affected only by the mutual friction (the 
corresponding characteristic equation is of order 6) . The CM modes, which start 
for Q — > at 1.03, are unaffected by the mutual friction. The modes which 
start for SI — > at 0, correspond to the rotational frequency of the spheroid 
in the low- J? limit. The right panel shows the same modes but when v/$(2 = 
4. Interestingly, while the viscous dissipation considerably affects the relative 
modes, its effect on the real frequencies of the CM modes is marginal. 

The damping via mutual friction, as seen from the left panel of Fig. 1, is 
maximal for rj/psuj = 1 and decreases to zero for ij/psto — * and rj/pgLU — > oo. 
This behavior is specific to the coupling between the superfluid and the normal 
fluid via the vortex state; the communication between these components is fastest 
when the magnitude of the forces on the vortex exerted by the superfluid and 
normal components are close. In the limiting cases the vortices are locked either 
in the superfluid (rj/psLu — > 0) or the normal fluid (r\j pguo — > oo) and hence the 
damping is ineffective. 

To conclude, the transverse-shear modes for the relative and CM modes re- 
main stable along the entire sequence of superfluid Maclaurin spheroids. 



4.4 Toroidal modes (I — 2, m — |2|) 

These modes correspond to \m\ — 2 and the motions in this case are confined to 
planes parallel to the equatorial plane. The toroidal modes are determined by 
the even in index 3 components of Eqs. ( |99| ) and (100) for the virials V-j, Vi-j, 
Ui-i and Ui-j, where i,j = 1,2. These equations can be manipulated to a set of 
four equations, which read 

(A 2 + 2f N uX + 4B 12 - 2Q 2 ) V 12 + flA(Vh - V 22 ) 

-2f N f s u\U 12 = 0, (119) 
(A 2 + 2f N uX + 4B 12 - 2Q 2 ) (V u - V 22 ) 

-AQ\V l2 - 2f N f s uX{U n - U 22 ) = 0, (120) 
(A 2 + 2J7/3A + 2/s^A + AA X - 2fi 2 \ U 12 



+0(1 - p')\{U n - U 22 ) - 2v\V l2 = 0, (121) 
"A 2 + 2i7/3A + 2f s v\ + 4Ai - 2Q 2 ) (U n ~ U 22 ) 



-4/2(1 - P')\U 12 - 2v\{Vu - V 22 ) = 0. (122) 

In the inviscid limit equations above decouple into separate sets for the CM and 
relative oscillations. The CM oscillations are described by the equations 

(A 2 + 4B 12 - 2Q 2 ) Via + OX{V xl - V 22 ) = 0, (123) 
(A 2 + 4B 12 - 2f2 2 ) (Vu -V 22 )- 4n\V 12 = 0, (124) 



The tensor virial method 25 








Fig. 2. The CM and relative toroidal modes of superfluid Maclaurin spheroid. The 
imaginary parts shown in grey are magnified by factor 10. Other conventions are the 
same as in Fig. 1. Note the secular instability at the bifurcation point e = 0.81 and the 
dynamical instability at the point e = 0.93. 



and their solutions are documented in EFE. The relative oscillations are de- 
scribed by the following equations 



(A 2 + 2[2f3\ + 4Ai - 2i7 2 ) (U n - U 22 ) - 4rtA(l - j3')U 12 = 0, (125) 
(A 2 + 2f2(3\ + 4Ai - 2^ 2 ) U 12 + flA(l - (3'){U n - U 22 ) = 0. (126) 
and the characteristic equation for the relative toroidal modes is: 



A 4 + 4/3/2A 3 + (8Ai + 4/3 2 l2 2 - 8/5'r? 2 + 4/3' 2 J? 2 )A 2 

+8/3i7(2A 2 - f2 2 )\ + 4(2Ai - f} 2 ) 2 = 0. 

In the frictionless limit the modes can be found analytically from 



(127) 



(A 2 + AAi - 2i7 2 ) 2 + AQ 2 \ 2 



(128) 
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which is factorized by writing A = ia. The two solutions are then 

c-1,2 = Q± \f±A\ - fi 2 - (129) 

and there are two complementary modes which are found by substituting — Q 
for fl. If the mutual friction is included the characteristic equation describing 
the relative oscillation modes is of order 4; in the presence of viscosity of the 
normal component, again the CM and relative oscillation modes couple and 
the characteristic equation is of order 8. The real and imaginary parts of the 
dissipative toroidal modes are shown in Fig. 2, for the same values of parameters 
as in Fig. 1. As for the transverse-shear modes, the CM and relative modes start 
at 1.03 and 1.63, respectively, when J? — > 0. In the inviscid limit (left panel), the 
CM modes are unaffected, as they should, while the relative modes are driven 
against each other and merge in the limit of strong coupling. The damping of 
the relative modes is finite, while for the CM modes it vanishes exactly up to 
the point of the onset of the dynamical instability at e = 0.95; beyond this 
point a mode becomes dynamically (i.e. in the absence of dissipation) unstable. 
If the kinematic viscosity is finite (right panel), the real parts of the relative 
modes are strongly affected, while the effect of the viscosity on the CM modes is 
marginal. The imaginary part of a CM mode changes its sign at the bifurcation 
point, where 2B 12 — ^ 2 and e = 0.813. This signals the onset of the classical 
secular instability of the Maclaurin spheroid. The new feature here is that the 
mutual friction contributes to the secular instability of a CM-mode. On the 
other hand, ordinary viscosity does not drive the relative modes unstable, in 
agreement with the fact that there are no neutral points for these modes along 
the entire Maclaurin sequence. One may conclude that the agents which break 
the superfluid/normal fluid symmetry can not cause an instability of the relative 
modes. The only possibility that the relative modes become unstable is a shift 
of the balance between the kinetic and potential energy perturbations, which 
might occur for compressible fluids, e.g. polytrops. This problem will be studied 
elsewhere. 



4.5 Pulsation mode (I — 2, m — 0) 



Pulsation modes (or breathing modes) are the generalization of the radial pul- 
sation modes of a sphere to the case of rotation. They correspond to I = 2 and 
m = indexes in the expansion in spherical harmonics. The pulsation modes 
arc determined by the full set of equations which are even in index 3. By suit- 
able combination of the equations for the virials Vi-i, Vi-j, Ui-i and Ui-j, where 
i = 1,2,3 and j = 1,2 the original set of equations can be reduced to 

(A 2 /2 + f N v\ + TBii - 2B 13 - Q 2 ) {V u + V 22 ) 



A 2 + 2f N v{l - 7 )A + 6B 33 - 2B 13 



^33 



+2Q\{V h2 - V 2 a) - fw/sfKUu + U 22 ) + 2f N f s v(l - j)U 33 
A 2 (^i ;2 - Vaji) - n\(v n + V 22 ) = 0, 



(130) 
(131) 
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A72 + 2A X - /2 2 + /<?j/A + Q/3X (Uu + U 22 ) 



U 



- A 2 + 2Q(3"X + 4A 3 + 2/51/(1 - 7 )A 

+212(1 - /3')A(C/i ;2 - £7 2;1 ) - uX{V xx + V 22 ) + 2i/(l - 7^33 = 0, (132) 

(A 2 + 2flj9A)(C7i.2 - U 2; i) - 0(1 - p')\(U u + U 22 ) = 0. (133) 

The characteristic equation is found by supplementing these equations by the 
divergence free constraint on the virials of the CM and relative motions: 



(134) 



i=l 1 t=l 1 

An equivalent form of the divergence free constraint for Maclaurin spheroids can 
be written in terms of eccentricity, (Vn + V 22 )(l — e 2 ) + V33 = and similarly 
for Un. 

In absence of viscosity the equations above decouple into two independent 
sets for CM and relative oscillations. The CM oscillations are described by 

(A 2 /2 + 4Bu - 2Bi3 - n 2 ) (V u + V 22 ) 

[\ 2 + 6B33 - 2B 13 ]v 33 + 2QX{V X ,2 - V 2 ,x) = 0, (135) 

(136) 



a 2 (Vi ;2 - v 2 , x ) - nx(v n + V 22 ) = 0, 



and coincide with the pulsation modes treated in EFE. The relative oscillation 
modes are defined by the equations 

A 2 /2 + QPX -Q 2 + 2A^j (U n + U 22 ) 

+2QX(l - p')(U ll2 - U 2 .i) - (A 2 + 4A 3 + 2f2~p"X)U 33 = 0, (137) 



A 2 + 2/2/3A) (U h2 - t/ 2; i) - nX(l - p')(U n + U 22 ) = 0, 



(138) 



which can be combined to: 



(A 2 + 2f2f3X - 2fl 2 + 4AA (A 2 + 2f2{3X) + 4f2 2 X 2 (l - (3') 2 (U n + U 22 ) 



[x 2 + 2J7/3A) [x 2 + 2J7A/3" + 4A-_ 



U 



33 



0. 



(139) 



The mode s ar e found by supplementing this equation by the divergence free 
condition ( 134 ). The third order characteristic equation in the inviscid limit is 

(3 - 2e 2 )A 3 + [8(3n + 4(3" Q - 40 + (3")e 2 n]X 2 + [4Ai + 8(1 - e 2 )A 3 
+2f2 2 + 4(3 2 f2 2 - 8/3' Q 2 + 4(3' 2 fl 2 + 8(1 - e 2 )(3(3" !7 2 }X 



+8A 1 (3H + 16(1 - e 2 )A 3 pn - 4/3J7 3 



0, 



(140) 



where the trivial mode A = is neglected. In the frictionless limit we find (A = ia 
as before) 

1/2 



a = ± 



2i7 2 + 4^i +8^3(1-6^ 
(3 - 2e 2 ) 



(141) 
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The pulsation modes for a sphere follow in the limit (e, fl) — > 0: for a sphere 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

e e 

Fig. 3. The CM and relative pulsation modes of superfluid Maclaurin spheroid. Con- 
ventions are the same as in Fig. 1 and 2. 

Ai/(irpG) = 2/3, and Eq. ( Jl4l| ) reduces to a 2 = 8/3 [a is given in units of 
(npG) 1 / 2 }. This result could be compared with the pulsation modes of an or- 
dinary sphere: a 2 — 16/15. Thus a superfluid sphere, apart form the ordinary 
pulsations, shows pulsations at frequencies roughly twice as large as the ordi- 
nary ones. In the general case where the viscosity of the normal fluid is taken 
into account the characteristic equation is of fifth order. The real and imaginary 
parts of the roots are shown in Fig. 3. In the inviscid limit the CM modes are 
again unaffected, while the relative modes are suppressed by the mutual friction. 
The damping of these modes is maximal when rj/ psoJ = 1 and the motions cor- 
respond to stable, damped oscillations. In the presence of viscosity, the relative 
modes are strongly damped and eventually become neutral. The CM modes are 
weakly affected. The imaginary parts remain always positive, i.e. one finds sta- 
ble, damped oscillations. Note that at the point where a relative pulsation mode 
becomes neutral the number of the imaginary components increases by one (as 
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generally expected for the roots of polynomials), but there are no changes in the 
sign of the imaginary part of the neutral mode. 

5 Summary 

Let us briefly summarize the main qualitative features of the oscillations of 
superfluid self-gravitating systems : 

• The oscillation modes of the superfluid ellipsoids separate into two generic 
classes which correspond to co-moving and relative oscillations. The oscil- 
lation frequencies of these two classes have distinct values in the both slow 
and rapid rotation limits. The first class of oscillations is identical to those 
of classical single-fluid ellipsoids in the incompressible and inviscid limits. 
Corresponding modes are undamped if the Euler equations of fluids are 
symmetric/anti-symmetric with respect to the interchange a <-> f3. When 
the fluids are coupled by mutual friction and mutual gravitational attrac- 
tion this symmetry is preserved. 

• The second class of oscillations, which is new, corresponds to relative mo- 
tions of the fluids. The modes are damped by the mutual friction between 
the superfluid and the normal fluid. These modes correspond to stable oscil- 
lations. 

• The co-moving (CM) modes emit gravitational radiation and undergo radi- 
ation reaction instabilities in full analogy to single-fluid ellipsoids |13j. The 
relative modes do not emit gravitational radiation at all, since the mass cur- 
rent associated with them is zero. This picture must hold true for a more 
general class of Chandrasekhar-Friedman-Schutz (CFS) radiation reaction 



instabilities, which are intrinsic to self-gravitating Newtonian fluids 14 . 
• If the a «-> (3 symmetry is broken the two classes of modes mix, for ex- 
ample, when the normal fluid is viscous. The main effect of the mixing is 
the renormalization of the mutual friction and viscosity. The relative modes 
remain stable as there are no distinct neutral points for these modes along 
the ellipsoidal sequences. The CM modes become unstable at the classical 
points of onset of secular/dynamical instabilities, for example at the point 
of bifurcation of the Maclaurin spheroid into a Jacobi ellipsoid. 

These qualitative features are based on general symmetries of underlying hy- 
drodynamic equations and the conditions of equilibrium of self-gravitating fluids 
which are independent of the superfluid nature of underlying fluids (e.g. the 
existence of bifurcation points) . Therefore we may conclude that these features 
will be preserved in more complex models of oscillations of superfluid neutron 
stars. 
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